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Magnetic subbands of graphite have been investigated by magneto-infrared reflectance spec- 
troscopy at 4K in fields up to 31T. Both Schrodinger-like (K-point) and Dirac-like (H-point) Landau 
level transitions have been observed, and their magnetic field dispersion are analyzed by a newly- 
derived limiting case of the Slonczewski- Weiss-McClure model. The values of the band parameters 
are evaluated without using sophisticated conductivity peak lineshape analysis. In this work, several 
less-explored band parameters are determined from the experimental results and they are known to 
result in electron-hole asymmetry and the opening of an energy gap between the electron and hole 
bands in multilayer and bilayer graphene systems. 



The tight-binding model of graphite's band structure, 
first calculated by Wallace [l| and later extended by Slon- 
czewski, Weiss, and McClurejl], has enjoyed a renais- 
sance in recent years due to its successful application in 
understanding many properties of single and few-layer 
graphene. 0, Hf The model, known as Slonczewski- Weiss- 
McClure (SWMc) band theory, rests on seven tight bind- 
ing parameters, 70 to 75 and A, characterizing interac- 
tions between near neighbor atoms in the lattice. The 
physical representations and estimated values of these pa- 
rameters have been illustrated and organized by Dressel- 
haus et al.\$ and Zhang et a/.@|. While the SWMc model 
has been used in understanding graphitic materials for 
the last 60 years the values of less-explored band 

parameters (74, 75 and A) remain uncertain. As these 
parameters may play an expanded role in interpreting ob- 
served asymmetries in the electronics of graphene, there 
is an increased incentive to accurately measure them. 

Infrared (IR) (magneto)-spectroscopy is a powerful 
tool in resolving energies in the band structure of 
graphitic materials, and has already been employed to 
study the unique electronic states of massless Dirac-like 
fermion in monolayer graphenell* 17 1, massive chiral 
fermions in bilayer graphene 7l | 22j -24|, and rich spec- 
tral features in graphitejfjl [l2 7fl3 7^18 - 2 1 1 . where effects 
attributed to both massless holes and massive electrons 
have been seen. The results in graphite have been mod- 
eled by assuming that its band structure is a combination 
of a monolayer model, where the holes at the Brillouin 
zone H-point behave like massless fermions, and an ef- 
fective bilayer model (EBM) [191-El] . where the electrons 
at the K-point behave like massive fermions with an ad- 
justed coupling constant between layers. However, the 
EBM ignores the band parameters governing the inter- 
layer hoppings between unstacked and stacked/unstacked 
sublattices (73, 74), next-nearest interlayer hoppings (72, 
75), and the on-site energy difference between the 
stacked and unstacked sublattices (A). Excluding these 



parameters removes the possibility of electron- hole (e— h) 
asymmetry, which has been observed to be significant in 
bilayer graphene and graphite. 

The interlayer hoppings can also drastically change the 
band structure near the charge neutrality point (CNP), 
i.e. the point where the electron and hole bands touch, 
and lift the degeneracy of the doubly degenerate e — h 
mixed Landau level, (LL-1) and (LLO). Here, we use the 
level indexing notation of ref. [ll[. The (LL-1) and (LLO) 
refer to the two LLs near the CNP, and the term, hn— 
em, refers to a transition from the n-th hole LL to the ru- 
th electron LL. The splitting of the two e — h mixed LLs 
was observed recently by Li et a/. [18], and the transition 
between these two states exhibits a potential 3D to ID 
crossover due to graphite's anisotropic properties. 

In this communication, we use IR spectroscopy to mea- 
sure the magnetic-field dispersion of graphite's LLs and 
fit our measurements to a newly-derived limit of the 
SWMc model that includes higher-order band parame- 
ters. The intense magnetic fields applied distinguish the 
two e—h mixed LLs and split of the e — h symmetry of 
several interband transitions. These lifted degeneracies 
are accounted for by our modified SWMc model, and val- 
ues for the set of band parameters are found consistent 
with those previously reported. @, 0] 

Our measurement set-up uses Kish graphite flakes sta- 
bilized on the Scotch tape and placed in a cryostat held 
at 4K subjected to magnetic fields of up to 31T. IR re- 
flectance spectra are measured by a Fourier transform IR 
interferometer using light pipe optics. The reflectance is 
then normalized by taking the ratio of the sample spec- 
tra to that of a gold mirror. The acquired reflectance 
spectra are modeled by Drude-Lorentz functions using 
RefFit. Sharp or minute features attributable to noise 
were smoothed prior to modeling. To show the quality of 
the modeling, the smoothed reflectance spectra are plot- 
ted in Fig. 1 overlaid with the original spectra at every 
other Tesla. Further analysis reveals that the observed 
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FIG. 1: Reflectance spectra of graphite: The spectra are ver- 
tically shifted for clarity. The smoothed (original) spectra 
are shown in thick red(thin green) lines. Two dashed lines are 
used as guides for spectral features following \AB-dependence. 



oscillatory spectral features are attached to the K- and 
H-point transitions, exhibiting linear in B and \B- de- 
pendence, respectively. The features resulting from e — h 
mixed LLs are located at low frequencies, while e — h 
asymmetry results in splitting in several of the interband 
transitions at high frequencies. 

The dynamic conductivity, shown in Fig 2, can be cal- 
culated from the modeling of Fig. 1, and exhibits a se- 
ries of magnetic-field dependent conductivity peaks. The 
overall profile of the dynamic conductivity is consistent 
with a recent theoretical calculation in the presence of 
a moderate magnetic field pit], and most K-point transi- 
tions show sizeable deviation from the linear dependence. 
The two lowest-energy modes, resulting from the (LL-1)- 
(LLO) and (LLO)-el transitions, are displayed in Fig. 2 
(a). It is obvious that magnetic- field dispersion of the 
(LLO)-el transition is not linear, instead exhibiting an 
interesting crossover from linear (low-field) to nonlinear 
(high-field) dependence. 

Graphite has been shown to exhibit a universal con- 
ductance Go = e 2 /Ah for photon energies between 0.1 
to 0.6eV at zero magnetic field. 26s] (leV= 8065cm -1 ) 
The sheet conductance of our sample, calculated using 
G(ui) = co<j(ui) with interlayer distance en = 3.35A, is 
shown in the inset of Fig. 2 (b). The calculated value is 
found larger than Go = e 2 /Ah, but closer to Ae 2 /h. Ac- 
counting for this discrepancy, one should note that the 
conductivity is not calculated from a true Kramer-Kronig 
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FIG. 2: Real part (m) of the dynamic conductivity: The 
curves are shifted vertically for clarity. For comparison, o\ at 
B — 0T is shown in a dotted line. To help the readers locating 
the modes in (b), two vertical lines mark the peaks associated 
with the (LLO)-el and h3-e2 transitions. Two guide lines 
similar to Fig. 1 follow the H-point transitions. Inset: Sheet 
conductance per graphene layer at B = 0T and 31T are shown 



in units of Go = e /Ah. 



transformation, but a result of modeling the reflectance 
spectra using finite numbers of Drude-Lorentz functions. 
As one cannot precisely assign Lorentz modes outside 
the measurement range, the optical weight for the modes 
located within the measurement range tends to be over- 
estimated. Moreover, the incident angle and the mis- 
alignment (between the sample and the gold reference) 
can also affect the result. On the other hand, the uni- 
versal conductance is dependent on the short-range dis- 
order, "ripples" in graphene's atomic structure and the 
presence of charged impurities 27], |28j, all of which are 
likely to be present in Kish graphite. The focus of this 
work is the magnetic-field dispersion of the LLs, which 
are not impacted by the frequency-independent universal 
conductance. 

Mode energies extracted from the conductivity peaks 
of Fig. 2 are plotted in Fig. 3, with modes following \f~B- 
dependence shown in (a) and those following linear de- 
pendence shown in (b). The (LLO)-el transition follows 
linear dependence below 12T, but \f~B- dependence at 
higher fields, and thus is plotted in both panels. Anoma- 
lies in mode energies can be found around the charac- 
teristic phonon energies of gra phitejlillO], particularly 
along several LL transitions (H: -1 to 0, K: (LL0)-e2, 
hl-e2), which will be further investigated elsewhere. 

In order to invoke the SWMc model, we start with the 
theoretical observation that the LL structures for 73 = 
and 73 ^ are similar at high fields III, |2l| , and so set 
73 = 0, which ignores the trigonal warping in the band 
structure. The energy of the magnetic subbands can then 
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be obtained using the following self-consistent Eq.jgj], 
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■gp^r-], in which cq and ao are lattice parameters. In 
SWMc model, Ei, E 2 , and E 3 can be represented asQ, 



E h2 = A ± 7l r + i 75 r 2 and e 3 = i 72 r 2 . 

At H-point, where r = 0, the levels areQ, 
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FIG. 3: Transition energies as a function of magnetic field: 
Three green dotted horizontal lines mark the characteristic 
phonon energies in graphite. Black solid lines are fittings 
using Eq. (3) for panel (a) and Eq. (6) for panel (b). The 
purpose of the dashed lines is explained in the text. Inset: 
Schematic energy diagram at the H- and K-point. 



Using the modes following V B-dcpendence, we find that 
the data agree well with the Fermi velocity c = 1.03 x 10 6 
m/s (70 = 3.18cV) and assign the two highest energy 
modes to the n = — 1 to n = and n = — 1(— 2) to n = 
2(1) transitions. To see if the parameter A is finite, the 
energy E—± t a for the n = — 1 to n = transition is scaled 
by a factor of 1 + ^/2, and plotted in a dotted line in Fig. 3 
(a). If A is finite, the scaled energy (1 + V2)i?_i,o would 
be larger than the interband transition energy E-i 9.(l5l] 
From the difference, |A| is estimated to be 60 cm -1 . 

For K-point transitions, the levels follows a linear in B 
dependence when E is small, and can be written asQ 



[±\jn(n + 1)lo 2 + -5uj 2 



(n + -)8w]qB, (4) 



in which u) determines the energy separation between the 
LLs and Suj results in the e — h asymmetry. It is clear 
from the data that most K-point modes are not linear 
at high fields. Though having a nonlinear magnetic field 
dispersion, the EBM[2(J [2lj ignores the e — ft asymme- 
try; an effect one cannot sufficiently account for by using 
two different Fermi velocities [21], since the correction is 
apparently dependent on the LL index n. On the other 
hand, It worths noting that the EBM[20|, HU and the 
bilayer model[3l| (r = 1) can be obtained by approxi- 
mating the self-consistent Eq. (1) of the SWMc model 



as, 
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and use a trial e 2 = 2 " 7 I ~ 1 E 2 on the right-hand side. 
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To derive a working formula to account for the non- 
linear magnetic field dependence of K-point transitions 
at high fields and with e — ft asymmetry, it is instruc- 
tive to re-examine the Eq. (1). If the term (-^-) 2 is 
neglected, Eq. (1) becomes a simple quadratic equation. 

Using graphite's band parameters, the term (-^§-) 2 can 
be neglected when -g >> 8 cm _1 /T, a condition can be 
met for any n > 1 levels. With this approximation, the 
magnetic levels can then be calculated to be 
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Interestingly, when v = 0, Eq. (6) can be written as 
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: (n+1/2), and uj c is the cy- 
clotron frequency. The first term in Eq. (7) agrees with a 
solution of the Dirac equation in a magnetic field except 
for that the value of n' was supposed to be integers. [32j 
Eq. (7) mimics the relativistic kinetic energy of a mas- 
sive particle. Using this analogy, the effect of v may 
be regarded as a Doppler effect, which causes blueshifts 
(rcdshifts) to electrons (holes). 

In analyzing the e — ft asymmetry, potential assign- 
ments for the higher-energy modes can be arbitrary (due 
to the large number of possible combinations), thus we 
use only the lower-energy modes to determine the band 
parameters. Good agreement with the data can be 
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achieved using uiq = 14.5 cm~ 1 /T, Sujq = —3.4 cm _1 /T 
and v = 0.05. The calculated transition energies are 
shown in solid lines in Fig. 3 (b) and the higher-energy 
modes are assigned to the transitions that best describe 
the magnetic-field dispersion of the transition energies. 
Near the K-point, i.e. jiT » A + ^ 5 T 2 - , the 

fitting parameters, uj and buj, are related to the band 
parameters as 



1 

7iT' 
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(8) 



From the values of fitting parameters, we extract 71 = 
0.38eV, 74 = 0.08eV and -72 + 75 = 0.048eV. It is gen- 
erally agreedjl, 0] that 72 = — 0.02eV, so 75 is estimated 
to be 0.028eV. 

For the (LL-1) and (LLO) levels (i.e. n = 0), one can 
solve the cubic Eq. in Ref. [l(| to find 



e = and e ~ —SujqB, 



(9) 



respectively. Since the (LLO) level will be partially oc- 
cupied by electrons, the transition between the (LL- 
1) and (LLO) levels are made by those electrons away 
from the K-point. The parameter 5oj(k z ~ kp)q is 
found to be — 3.7cm /T, which is slightly larger than 
5u>{k z = 0)q = — 3.4cm /T; therefore, the sublatticc 
energy difference A is likely to be positive. 

The (LLO)-el transition exhibits an unusual crossover. 
At lower fields, its energy agrees with the prediction 
of the linear SWMc model with a slope of 17cm _1 /T. 
(dashed line in Fig. 3 (b)) At higher fields, the peak 
broadens quickly with increasing magnetic field, and the 
magnetic-field dispersion becomes Dirac-like with a slope 
of 54cm~ 1 /T 1/2 - (dashed line in Fig. 3 (a)) A simi- 
lar behavior has been reported previously [18j. where the 
(LLO)-el mode has a larger slope for B < 8.5T, but a 
smaller slope for B > 8.5T. In addition, we find that 



(ll-i)-(lw) + E (LL0)-ei < E (LL-i)-ei, a discrepancy 



that cannot be accounted for by simple pictures assum- 
ing LLs are independent with each other. Except for 
the (LLO)-el transition, all other transitions are well de- 
scribed by the modified SWMc model (Eq. (6)), so the 
underlying mechanism for this anomaly is specific to the 
(LLO) level. Anomalous phenomena have been observed 
pertaining to this level, and it has been conjectured that 
a charge density wave [33H351 has an effect specific to 
it, though it should be noted that these anomalies were 
found at lower temperatures and at higher fields. 

In summary, magnetic subbands in Kish graphite 
are investigated using magneto-infrared reflectance spec- 
troscopy up to 3 IT. The high magnetic field distinguishes 
the two e — h mixing LLs near the CNP and the splitting 
of the interband transitions. We propose a new limit of 
the SWMc model to describe the K-point transitions and 



use it to estimate the values of the band parameters, 

70 = 3.18eV, 71 = 0.38eV, 72 = -0.02eV, 73 = 0, 
74 = 0.08eV, 75 = 0.028eV, and |A| = 0.008eV. 
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